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Detailed derivation of the master equation and the corresponding time evolution of the cavity 
radiation of a coherent beat laser when the atoms are initially prepared in a partial coherent super- 
position is presented. It turns out that the quantum features and intensity of the cavity radiation are 
considerably modified by the phase fluctuation arising due to the practical incapability of preparing 
atoms in the intended coherent superposition. New terms having an opposite sign with the contri- 
bution of the driving radiation emerged in the master equation. This can be taken as an indication 
for a competing effect between the two in the manifestation of the nonclassical features. This, on 
the other hand, entails that there is a chance for regaining the quantum properties that might have 
lost due to faulty preparation by engineering the driving mechanism and vice versa. In light of this, 
quite remarkably, the cavity radiation is shown to exhibit nonclassical features including two-mode 
squeezing and entanglement when there is no driving and if the atoms are initially prepared in 
a partial maximum atomic coherence superposition, contrary to earlier predictions for the case of 
perfect coherence. 

PACS numbers: 42.50.-p, 42.50.Ar, 42.50.Gy 



I. INTRODUCTION 



In a two-photon three-level laser, the nonclassical fea- 
tures of the radiation are predominantly attributed to the 
atomic coherence. It is evident that the atomic coherent 
superposition in the cascade scheme can be induced via 
coupling the levels between which a direct transition is 
electric dipole forbidden, apparently the upper energy 
level | a) and the lower energy level |c), by external radia- 
tion with high amplitude PHD] . The same can be achieved 
by pumping the atoms with external magnetic field when 
magnetic dipole transition is allowed. The coherence can 
also be induced by preparing the atoms initially in ar- 
bitrary coherent superposition of these levels & 14 or 



by using these mechanisms simultaneously p^TfT6|. In 
the nondegenerate configuration, when the atoms spon- 
taneously decay from \a) to |c) via the intermediate en- 
ergy level \b), two photons with different frequencies are 
generated. For the sake of clarity, the amplification of 
light when spontaneously emitted photons are correlated 
due to atomic coherence induced via the initial prepa- 
ration and external driving mechanism is designated as 
coherent beat laser. In order to exhibit the required las- 
ing, the atoms are assumed to be injected into a doubly 
resonant cavity at constant rate and readily left after 
spontaneously decayed to levels that are not involved in 
the process. 

Despite frequent assumption, current practical capa- 
bility is not in a position to yield arbitrary coherence 
as expected (the review of the arbitrary atomic coherent 
preparation is found for instance in [ItJ )• Essentially, the 
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laser to be deployed in the preparation has a bandwidth 
(along with dissipation processes) that leads to fluctua- 
tions in the prepared coherent superposition. Quite ob- 
viously, there is also a time lag between the preparation 
and injection of atoms into the cavity let alone the in- 
consistency in injection mechanism that causes change 
in the atomic coherence. Hence the practical limita- 
tion of preparing the atoms in the intended arbitrary 
coherent superposition acclaims the issue of the partial 
prepartion of the same as an important aspect of the 
system under consideration. In this line, recently Qamar 
et al. 8] have studied three-level cascade laser with in- 
jected phase fluctuation and predicted that it influences 
the entanglement generation strongly. With this back- 
ground, the effects of the phase fluctuations on the in- 
tensity and nonclassical features of the cavity radiation 
when the atomic coherence is partially prepared initially 
and then induced by external driving mechanism latter 
(coherence beating) would be investigated in this paper. 
In earlier communications, efforts have been made to in- 
corporate the phase fluctuations arising from the pump 
radiation as well (which is usually designated as pump 
phase) [J, H, [H, [l9| . For the sake of bringing the com- 
plication that otherwise would be surfaced to the man- 
ageable level, the contribution of the pump fluctuation 
(even though its effect is anticipated to be significant un- 
der some circumstance) is not taken into consideration. 
Moreover, the initially prepared atomic coherent super- 
position is presumed to decay at arbitrary rate due to 
physical phenomena like vacuum fluctuations |15l |17| . In 
one of recent contributions, the cavity evolution of the 
latter with arbitrary initial coherent superposition has 
been considered where significant enhancement in gener- 
ated entanglement for certain cases of interest has been 
reported [15|- It is hoped that bringing these competing 
processes together may provide a solid background to ac- 
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tually attest the potential of this system as a source of a 
bright robust entangled light. 

Furthermore, according to Caves classification [20| . 
two-photon three-level cascade laser falls under phase- 
sensitive amplifiers HI| • This can be evinced by re- 
sponding differently to the quadrature phases in the from 
of unequal gain or unequal noise or both. In view of this, 
the partial initial prepartion of the coherent superposi- 
tion is believed to lead to modifications in correlation 
that have been characterized as strong and taken to be 
a basis for the prediction of substantial degree of two- 
mode squeezing jl j ■ entanglement 0, H, HH, [23 , violation 
of Cauchy-Schwarz inequality and violation of Bell- 
CHSH inequality [23| in various schemes. It, therefore, 
appears natural to expect a substantial change in these 
quantum features in connection to the present consid- 
eration. Since the consequent atomic transitions would 
be altered it would not be difficult to envision that the 
mean photon number would also be perceivably affected, 
eventhough serious investigation of the effect of the phase 
fluctuations on the intensity of the generated radiation 
is either lackingor unobserved in the master equation in 
earlier reports [H, H, H, [H, ■ As a result, in the present 
contribution, the effect of the phase fluctuations on the 
intensity of the generated radiation is also studied. 

To achieve these goals, the master equation and ex- 
pressions that represent the evolution of the cavity radi- 
ation in terms of the c-number variables associated with 
the normal ordering arc derived following the approach 
outlined in |15| and references there in. Although the 
scope of this work is restricted to deriving the time evo- 
lution of the cavity variables, it is basically possible to 
make a comparative study of the statistical properties 
and quantum features of the generated cavity radiation 
employing the obtained results. Contrary to the current 
trend of solving the resulting coupled differential equa- 
tions applying the numerical calculations, the analytical 
approach is chosen despite the involved rigor for the sake 
of avoiding the approximations that otherwise required. 



II. MASTER EQUATION 

Interaction of pumped nondegenerate three-level cas- 
cade atom with a resonant two-mode cavity radiation can 
be described in the rotating-wave approximation and the 
interaction picture by the Hamiltonian of the form 

H =ig[a\a) (b\ - \b) (a\a 1 + b\b) (c\ - \c) (bp] 

+ i^[\c){a\-\a){c\}, (1) 

where O is a real-positive constant proportional to the 
amplitude of the driving radiation and g is a coupling 
constant chosen to be the same for both transitions, a 
and b are the annihilation operators that represent the 
two cavity modes. Although the pumping laser has a 
bandwidth that consequently leads to phase fluctuations, 



this issue is not considered here for limiting the involved 
mathematical rigor. However, due to emerging various 
quantum effects, it is assumed that the atoms can only 
be initially prepared in a partial coherent superposition 
of the upper and lower energy levels. In view of this, the 
initial state of the three-level atom is taken to be 

|*a(0)) =C a {0)\a)+C c (0)e^\c), (2) 

where C a (0) = (a\V A (0)) and C c {0)e iv = (c|*a(0)> are 
probability amplitudes for the atom to be initially in the 
upper and lower energy levels, and <p is an arbitrary phase 
between the two states. It is assumed that ip is randomly 
distributed about a fixed mean phase ipo. In which case, 
the phase can be defined as ip = ipo + Sip where Sip is 
small random fluctuations around ipo (whose effect can 
be adjusted at will by proper choice of the phase of the 
cavity radiation Q). Therefore, with no further ado, we 
set ipo = and take ip as fluctuations about po = 0. 
Fundamentally, it is the effect of these fluctuations that 
is going to be addressed in the present work. 

In line with Eq. @ , the corresponding initial reduced 
atomic density operator would be 

Pa(0) = pi° tt >|a)(a| +pW e -^|a)<c| 

+ p^e^\c)(a\+p^\c)(c\, (3) 

where = |C a (0)| 2 , = C a (0)C*(0), pg? = 

C c (0)C*(0), and p^J = |C C (0)| 2 . It is straightforward to 

see that pha and pi? are the probability for the atom to 
be initially in the upper and lower energy levels, whereas 

pi?e ±I(/3 represent the initial atomic coherence. This in- 
dicates that the introduced phase fluctuations do not 
affect the populations at the beginning. However, the 
density operator for the cavity radiation plus a single 
atom injected into the cavity at time t% is represented by 
PAn{t,ti) in which t — T < tj < t, where T is the time 
at which the atoms are removed from the cavity. Thus 
the density operator that describes all the atoms plus ra- 
diation in the cavity when the atoms are continuously 
injected into it at a constant rate r a can be expressed as 

PAR(t) = r a f p AR (t,t')dt'. (4) 

Jt-T 

Replacing the summation over randomly injected atoms 
to integration in a similar manner is known for quite a 
long time 0,E1- 

On the other hand, with the assumption that the atom- 
radiation density operator can be decorrelated into the 
atom and radiation parts at a time when the atoms are 
injected into the cavity and when they just left the cavity, 
it is proposed that 

-TLPAR{t) =r a [p A {ti)-p A {t-T)]p{t)-i[H, p AR (t)], 

(5) 
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where p A {t) = Pa(0)- Moreover, taking the trace over 
atomic variables and using the fact that 



-^_Pab{t) = g(apbb ~ Paaa + pac&) - —Pcb - labPab, 



(12) 



Tr A (p A (0)) = Tr A {p A (t - r)) = 1 
along with Eq. (JTJ) , one can readily obtain 

= g\Pab^ - Cltpab ~ tipbc + Pbcb ] 

at L 

+ ap ba - p ba a + bp cb - pcbb] , (6) 

in which p a p = (pt\f>AR\P) with a,j5 = a, b, c. 

Furthermore, it is not difficult to see on the basis of 
Eq. © that 



dt 



P~ a p(t) = r a (a\p A (0)\P)p - r a (a\p A (t - r)\P)p 



i(a\[H, p AR {t)]\ 



lappa? , 



(7) 



where the last term is added to account for various atomic 
decay processes including atomic spontaneous emission. 
Assuming the atoms to be removed from the cavity after 
they successfully decayed to energy levels other than the 
intermediate or the lower implies that 

(a\p A (t-T)\P) =0. 

Consequently, making use of Eqs. (fTJ), ([3]), Q, and not- 
ing also that |o) and \b) are atomic states whereas a and 
b are cavity operators that do not act on these states, 
one readily finds 

^/V*(*) = -lafiPaf) + r a [p [ °}5 aa 5 af) + p^e^SaaScp 

+ p^e i n ac S aP +pi°j6 ac 6 c p]p(t) 

— g[a> Pa/jSab + b^pbfjSac - apbpSaa ~ bp c f]5 a b 
+ PaaaSbf) + p a bb8 c fj - p a ba) 5 a {3 - p ac b' Sbfi] 

~ 2 [Pc^aa — PafiSca ~ PaaScfj + Pac^ap] , (8) 

from which follows 



^Paa(i) = r a p<£}p(t) + g{ap ba + pabtf) 



n 



(Pac + Pea) - 



(9) 



dt 



Pbb(t) = -gi^pab + Pb a a - bp cb - / o bc 6 t ) - TbPbb, 



(10) 



dt 



Pcc(t) = r a p c °Jp{t) - g(tfp bc + p cb b) 

+ ^(pac + Pea) ~ TcPcc, 



(11) 



dt 



d A , . 

J t Peb(t) 



Pac(t) = r a p^Je l(p p + g(ap bc - pabb) 



O 



(Pc 



laePa 



(13) 



n 



-g(p ca a ~ Pcc& + b^pbb) + -^Pab - lebpeb, 



(14) 



where Ej = 7^ and 7y(^y) with i, j — a.b,c stand for the 
atomic decay rate and the rate of dephasing (the rate at 
which the atomic coherent superposition decays), respec- 
tively. As one may clearly observe, the contribution of 
the phase fluctuations appeared only in Eq. (| 13)) . 

The solution of these coupled differential equations and 
their complex conjugate would reveal how each popula- 
tion and various correlations (atomic coherent superposi- 
tion among the three levels) evolve in time. Even though 
this is not the target of the present contribution, it is 
worth mentioning that solving these equations renders 
the response of the atom in regards to various assump- 
tions which would be an attractive issue by its own right. 

Basically, in the good cavity limit where the cavity 
damping rate k is much smaller than atomic decay rates 
(r, and Jij), it is possible to employ the adiabatic ap- 
proximation scheme (Tlj . Confining to linear analysis, 
which amounts to dropping the terms containing g in 
Eqs. ©, JTDJ}, (|TT|). ([15]). and then applying the adia- 
batic approximation scheme result 



faPaaPit) - ^Pac ~ ^aPaa = 0, 



Pbb = o, 



TaPbSKt) + ^Pac ~ ^ePec = 0, 



(15) 



(16) 



(17) 



r a p^e-^p(t) - -(p cc - p aa ) - lacPac = 0, (18) 

where we set p ac — p ca on the basis that once the atoms 
entered the cavity a random phase between the upper and 
lower energy levels can be ignored. It should be noted 
that the considered phase fluctuation is solely attributed 
to the quantum processes at the time of preparation and 
injection, but not to the cavity dynamics. Nevertheless 
it is incontestable that due to various physical processes 
there could be phase fluctuations in the ephemeral co- 
herent superposition of the atomic levels while the atom 
traverses in the cavity. Owing to the practical limitation 
of probing into this effect, it is rather ignored at least 
for the time being. In essence, the linear analysis is re- 
quired so that the resulting differential equations can be 
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analytically solvable. Since the quantum features in this 
system are associated with the correlation induced in the 
cascading phenomenon rather than the nonlinear process 
as in the other quantum optical systems, the lineariza- 
tion approach still holds while the nonclassical properties 
of the radiation is studied. 
Now setting 



P (0) 

raa 



(19) 



with — 1 < 77 < 1 , it is not difficult to verify that 

Jo) 



= -^f 2 and p^ac — ^2 ■ Moreover, upon setting 



Pec - 2 

r Q = T b = T c = T and 7 ab = 7 ac = 7 cfc = 7, it is possible 
to get with the aid of Eqs. ([IS]), 1(17)1. and (JUJ that 



Pact 



r a p 



2r( 7 r + n 2 ; 



[ 7 r(i - n) 



(20) 



r a p 



2r( 7 r + r> 2 ) 



[ 7 r(i + n) 



Pa 



r a p 



2( 7 r + n 2 



Ty/l — r} 2 e 



Or? 



(21) 



(22) 



with p = p(t). It is vividly seen that the phase fluctua- 
tions introduced during atomic preparation modifies the 
subsequent populations and coherence, although the pop- 
ulations are initially unaffected. Further scrutiny reveals 
that p aa , Pec, and p ac exhibit oscillatory nature entirely 
emanated from the phase fluctuations. This can be cited 
as the indication for the dependence of the intensity and 
nonclassical features of the radiation on the initially in- 
troduced phase fluctuations. 

Next, making use of Eqs. (0, OH), <H3>, (EE)), 

([2"2")l, and applying the adiabatic approximation scheme 
once again, one gets 



Pah 



gr a p 



7 2 (4 + C 2 )(l + CC) 
+77 (C'C - 2 X ) - (2C + Qy/T^e-* 



+St[C'(l + C'C) + 3r?C' - (2 - C'OVW^] 



(23) 



Pcb 



gr a p 



7 2 (4 + C 2 )(l + C'C) 



[a[C'(i + CO 



-3r?C' + (2 - COVl 3 ^ 

-St [2 (C' 2 +x)-??(C'C-2x) 



(24) 



where C = (' = Sl/T and x = 7/T- It is worth 

noting that like /5 ac , /5 a {, and /5 C (, also display oscillatory 
nature. 

Furthermore, upon employing Eqs. ([5]), (1231) . and ([2"3)l. 
the master equation can be put in the form 

= — [2apa) — a)ap — pa) a] 



~2B 
1 (AC. 



2 V B 
AD+ 



+ 



2B 
AD_ 

2B 

AE + 

2B 
AE_ 



bpw — a) pa — wbp + aa' p 



[2a' pa — paa' — aa'p] 

26 / 56 t -pb r b- b r bp 
aa) 

bpo* — a/ pa — pwb + paa/ 
a' pw — V s a) p + bpa — abp 
a' pw — pb'a) + bpa — pab 



2B 

ACQ. + CQ 
2B 



b'a} p — ply <v — abp + pab , (25) 



where 



2r a g 2 



B = (4 + C 2 )(1 + C'C), 



C ± =2C' 2 + 2 X ±r / (C'C-2x), 



D, 



(2C' + C)v / T^7 2 e 



±2^ 



e± = 37?c'-(2-c'c)yr 



(26) 



(27) 



(28) 



(29) 



(30) 



The contribution of the cavity damping which corre- 
sponds to the coupling of the cavity modes with envi- 
ronment modes via the coupler mirror is incorporated 
following the usual standard approach (2(| H?} • It is not 
difficult to observe that this master equation has a fun- 
damental difference in the form from earlier reports of 
various similar schemes [l5| . 

Based on the form of the master equation, it has been 
a customary designated for instance C+ as the gain of 
mode a and C_ as the lose of mode b. Even though the 
terms associated with C± have the same form as earlier 
results, in Eq. (f25)) . some of these terms are unambiguous 
included in D± as well. As a result, it is not appropriate 
concluding that C± is the gain or lose outright. Rather 
such jumbling implies that the mean photon numbers of 
the respective modes come from the contribution of terms 
related to C± and D± that has to be thoroughly ad- 
dressed. Since D± depend on the phase fluctuations, this 
result asserts that the intensity of the radiation should 
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depend on the same. This can be related to the fact that 
the phase fluctuations initially introduced unequivocally 
alter the subsequent atomic transitions responsible for 
the generation of the photons. With all limitations, the 
change in the transient populations is believed to be the 
main cause for the dependence of the mean photon num- 
ber on the initial preparation. In this line, although nu- 
merical calculations indicate that the mean photon num- 
ber depends on the phase, this fact has not been directly 
reflected in the master equation 0, H, [j| [TH, . More- 
over, it can readily be seen from Eq. (|2"51) that D± = 0, 
if f2 = or r\ = 1, in earlier reports unfortunately cither 
of these cases has been taken which might be the source 
of the disparity. 

On account of recent discussion that the degree of two- 
mode squeezing and entanglement decrease with the rate 
of dephasing [28], it can be observed that the statisti- 
cal properties and entanglement of the cavity radiation 
significantly depend on the rate at which the coherent 
superposition is decaying in this case as well. It is also 
noticeable that the last term in the master equation is 
solely related to the external driving radiation which we 
do not wish to dwell on it further here. Nonetheless, it 
may worth mentioning that the last term in the mas- 
ter equation appears due to the emergence of new terms 
that should have been zero when the phase fluctuations 
are not taken into consideration. On the basis that it is 
the mixed terms that are responsible for observing non- 
classical features, the effect of the phase fluctuations on 
the quantum features of the radiation enters into play 
in two ways. The first is associated with this emerging 
terms where the other is through the dependence of E± 
on ip. The exponential dependence of the quantum fea- 
tures as in E± has been touched upon earlier in different 
context. However, the contribution of the last term in 
the master equation is yet to be explored. To this effect, 
critical scrutiny reveals that the sign infront of wa^p and 
pab in E± and the last term of the master equation are 
different. This suggests that the external driving and ini- 
tial preparing mechanisms may have a competing effect 
in which the quantum features that would be lost due 
to incapability in preparing the atoms be compensated 
by the external driving mechanism. By and large, this 
would be an encouraging result if it is found to be correct. 



III. EQUATIONS OF EVOLUTION 

Quite often, the c-number Langevin equations are 
found to be easier mathematically to handle than the cor- 
responding operator equations. This is one of the reasons 
for choosing the stochastic differential equations over the 
operator equations that can be derived from the master 
equation directly. Moreover, recent works make it evi- 
dent that the stochastic differential equations associated 
with the normal ordering of the cavity mode variables are 
important tool in comfortably studying the quantum fea- 
tures of the radiation [3. l2l|. Therefore, in this section, 



these equations would be obtained applying the perti- 
nent master equation. To this end, employing Eq. (j2"5)l 

and the fact that gj(6(i)) = Tr (j§6\ (where 6 is any 

operator) it can be verified that 



dt 



[a(t)) 



2B 



[Bk + A(D_-C + )] (a(t)) 



A[ C{1 + (C)+ E_](b\t)), (31) 



dt 



(&(*)) 



2B 



[Bk + A(C- +£>+)] (&(*)> 



.4 



+ ^[C'(i + CC')-£+]« 



r (*)>> 



(32) 



dt 



(ot(t)o(i)) 



2B 



[2kB + A(D- +D+- 2C+)] 



x (at(t)a(i)) 

+ ^[E + +C(l + CO](a(t)b(t)) 



2^" 
A 



C'(i + c'C)]<aW(*)> 



+ —(2C + -D + -£>_), 



(33) 



x (b\t)b(t)) 

-2^ 



C(l + CC)](a(t)b(t)) 



-(a(t)b(t)) 



— [E + -C(l + COW(t)bHt)), (34) 



— [2kB + A(C_ - C+ + D_ + D+)] 



x (a(t)b(t)) 

-^[E + -C(l + CQ]{aHt)a(t)) 



2^" 
.4 



2B 



C'(1 + C'0]<& + (*)&(*)> 



(E+ - C'(i + C'0> 



(35) 



where a bar is put over D± and E± to indicate that 
stochastic average is taken over the phase since it is as- 
sumed to fluctuate about the central mean value (set to 
for convenience). From practical point of view, ad- 
dressing the contribution of every phase change seems to 
be unrealistic. In connection to this, assuming that the 
phase undergoes Gaussian random process [g, [Io|, fl8| . 
it would be more appropriate if the phase fluctuation is 
used instead of the actual phase. With this understand- 
ing and employing the fact that for Gaussian variables 

(exp±iSip) = exp — (5ip 2 /2) 
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exp±iip is replaced by exp— in the above equations. 
It is not difficult to realize that for Gaussian random 
process (6<p) would be zero, hence 8 — (S 2 (f>/2) represents 
the deviation which is generally designated as phase fluc- 
tuation. Basically the dependence of the properties of 
the cavity radiation on this parameter is the subject of 
present report. But it is worth noting that if it is possible 
to phase lock the laser employed in the initial prepara- 
tion and the fluctuations resulted due to other quantum 
processes can be negligibly small, taking the average over 
the fluctuations may not be strictly required. 

In the same manner, any other required time evolu- 
tion of the correlation can be generated using this mas- 
ter equation. It goes without saying that these coupled 
differential equations can be numerically solved, whereby 
the required information can be analyzed and the system 
under consideration can be studied. However, solving the 
same problem employing analytic approach is expected 
to provide a far more complete picture than the approx- 
imated numerical solutions can offer. In order to pave a 
way for solving this problem analytically, the operators in 
the above equations are put in the normal order. Hence 
the corresponding expressions in terms of c-number vari- 
ables associated with normal ordering for the first two 
equations (Eqs. ([3T]) and ([32]) ) can be rewritten in the 
form 

j t a(t) = -a + a(t) - b + (3*(t) + f a (t), (36) 
j t P{t) = -a-0(t) - b-a*(t) + f b (t), (37) 

where 

a± = \ + 4b [(2C ' + c)e T 2(c ' 2 + x) _ ^ C ' C _ 2x)] ' 

(38) 

b± = -4s [C'(l + CO ± [377C' ~ (2 - CC')©]] , (39) 



with = e~ fl yl — rf . No doubt that the effects of the 
phase fluctuation prominently enter via these coupled dif- 
ferential equations which believed to represent the dy- 
namical evolution of the system including the contribu- 
tion of the noise emanating from various quantum pro- 
cesses. With this conviction, the properties of the gen- 
erated radiation can be evaluated and then analyzed by 
directly solving these equations following straight alge- 
bra, even though the steps are somewhat lengthy. Even 
without going into details of the solution, it is clear to 
see that the actual effect of the phase fluctuation is as- 
sociated with the form of the parameters defined in Eqs. 

and (f3"51) . In this respect, critical scrutiny of these 
equations shows that the effect of the phase fluctuation 
appears in b± even when O = (no external pumping 
radiation). This indicates that the phase fluctuation re- 
mains to affect the mean photon number and nonclassical 



features of the radiation whether there is external driving 
mechanism or not, although its effect felt rather lower in 
certain cases when there is no driving mechanism. 

Moreover, f a {t) and /&(t) are the noise forces the prop- 
erties of which remain to be determined. For instance, 
the expectation values of Eqs. ([36j) and (|37|) would be 
identical to c-number equations corresponding to Eqs. 
(3D and (02D provided that (f a (t)) = and (f b (t)) = 
which implies that the noise forces have stochastic na- 
ture. The correlation between these noise forces can 
be evaluated by following somewhat tricky procedure in 
which the various combinations that emerge by apply- 
ing Eqs. (f36f and (|37|) are compared with the evolu- 
tion equations generated from the master equation where 
in the intermidiate step the formal solution of the same 
equations are taken. In the final evaluation of the cor- 
relations, the system is assumed to be unaffected by the 
noise force evaluated at latter time. Going through the 
outlined procedure reveals that 

(Ut')f a (t)) = I [2C 2 + 2 X + V(CC - 2 X ) 

-(2(' + 0®]S(t-t'), (40) 

</ 6 (i')/ 6 *(*)}=0, (41) 
+(2-('0®]5(t-t'), (42) 

(ti(t')fa(t)) = (fa(t')fa(t)) = (A(f)A(*)> = 0. (43) 

It is evident from these results that the effect of the noise 
source appears as it does in the corresponding quantum 
description. The properties of the correlations of these 
noise forces are related to the normal operator ordering 
which directly associated with the vacuum fluctuations 
in the environment and cavity modes. Comparing Eqs. 
(|40| and (|4l"|) reveals that the correlation properties of 
the noise corresponding to mode a and mode b are dif- 
ferent unlike the other quantum optical systems. This 
disparity basically reflects the difference in the number 
of photons in mode a and mode b [U [22| . With mini- 
mum effort, it is also possible to see from Eqs. (|4l7|) and 
(|4"2|) that the correlations of the noise forces depends on 
the phase fluctuation. This generally indicates that the 
mean photon number and the nonclassical features of the 
radiation considerably depend on this fluctuation. It is 
by now a common knowledge that the nonclassical fea- 
tures are witnessed in cross correlation of the noise forces. 
With this background, it is possible to predict that the 
phase fluctuation keeps on affecting, let us say, the gen- 
erated entanglement as long as the atoms are injected 
into the cavity and they are not initially prepared in ei- 
ther the upper or the lower energy levels (or as long as 
there is initial atomic coherent superposition). Extend- 
ing this argument in the same line shows that the effect 



of the phase fluctuation would be maximum when r\ = 0, 
that means, when the atoms are initially prepared with 
equal probability to be either in the lower or upper en- 
ergy level. In connection to the absence of entanglement 
in this case when O = [HI ], it is expected that the ex- 
ternal phase fluctuation initiates nonclassical correlations 
via indirectly creating electron pathway. Unfortunately, 
the effect of the phase fluctuation on properties of the 
cavity radiation in connection to the change in the am- 
plitude of the external radiation and rate of dephasing is 
not evident as one wishes it to be. 

In order to go deeper, it is desirable solving the prob- 
lem analytically. In this effect, it is straightforward to see 
that Eqs. (|36|) and (|37]) form coupled differential equa- 
tions that can be readily solvable following the approach 
outlined in [15] . Following straight algebra, one can ver- 
ify that 

a(t) = F+(t)a(0) + G+(t)j9*(0) + H+(t) + I+(t), (44) 



(3(t) = F_ (i)/3(0) + G_(t)a*(0)+ H-(t)+I-(t), (45) 
where 

F±(t) = e- At [coshrf±psinhrf], (46) 



G±(t) = -q±e~ M sinhet, 



(47) 



H+(t) 



H-{t) 



■ A (*-*''[coBhe(t-t') 



■psinke(t-t')]f a (t')dt', 



-^-^[cosheit-t') 



■psmhe(t-t')]f b (t')dt', 



(48) 



(49) 



I+(t) = -q+ e-W-^sinheit-t')^ (t')dt', (50) 



/_(*) = -g_ / e-^-t'hmheit-t^f^t^dt', (51) 
with 



AZ 
~2B' 



(52) 



A = \ + 4b [(2C ' + c)e " r ' (c ' c ~ 2x)] ' (53) 



2[C' 2 + X] 



-qi + C'OT^c- (2-CC)e] 

Q± = (55) 



in which 



\/C' 2 (l + CO 2 + 4[C' 2 + X? - [3t?C - (2 



C'C)e] 2 . 

(56) 



(54) 



It perhaps worth mentioning that Eqs. (|44|) and (|45[) 
along with the associated parameters are used to cal- 
culate various quantities of interest. It is noticeable that 
these solutions are well behaved functions at steady state 
provided that A — e > 0. As a result, the case for which 
A = e is designated as a threshold condition. As criti- 
cal scrutiny of the expressions following from the master 
equation reveals, this mathematical condition is directly 
related to the uncertainty condition 

Despite the fact that the form of the master equation is 
considerably shifted from results reported in earlier com- 
munications due to the introduction of the phase fluctu- 
ation, the general form of the time dependence of a(t) 
and (3{t) remains to be unaltered. This asserts that the 
auto and cross correlations of the radiation modes can be 
calculated in the same way as earlier cases. Even then, 
it would be essential noting that the form and content 
of F±, G±, H±, and I± are significantly changed. More- 
over, it is not difficult to see from Eqs. (fS"2"j) . ([53)l . ([51)1. 
(fS"5")) . and (|5"6")) that the parameters (e, A, p, and q±) that 
determine the content of the provided solutions depend 
on the initially prepared coherence, the amplitude of the 
driving radiation, the rate at which the atoms sponta- 
neously decay, and the phase fluctuation. In addition to 
these, the cavity damping constant, the rate at which 
the atoms are injected into the cavity, and the coupling 
between the atom and the radiation significantly affect 
the photon statistics and the nonclassical features of the 
generated radiation. Hence, it may not be difficult to 
comprehend that providing the effects of phase fluctu- 
ation may not be simple and straightaway as the form 
of the solutions might suggest. However, it should be 
emphasized here that any required correlation can be de- 
termined and consequently in depth analysis can be made 
using these solutions. Ironically, the investigation in this 
direction is quite an involving task. As a result, for the 
sake of clarity, the detailed study of the quantum prop- 
erties of the cavity radiation is differed to subsequent 
communications. 



IV. CONCLUSION 

In this contribution, detailed derivation of the evolu- 
tion of the cavity radiation of two-photon coherent beat 
laser when the initial atomic coherent superposition is 
partially prepared is presented. Generally, the master 
equation is found to have additional terms resulting from 
the change of sign due to the phase fluctuation associated 
with partial preparation. Particularly, the mixed terms 
that are closely related with nonclassical features of the 
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radiation turn out to have opposite sign with correspond- 
ing contribution of the driving radiation. Without going 
into details, this disparity alone suggests that the exter- 
nal driving mechanism and the partial initial preparation 
offers a competing effect. This, on the other hand, im- 
plies that there is a possibility for regaining the quan- 
tum properties that might have lost due to incapability 
of preparing atoms in a particular atomic coherence by 
externally driving them with classical radiation. Quite 
remarkably, this seems to be an obvious result when one 
reviews how each process affects the atoms. In connec- 
tion to this, critical observation of Eq. (pfHj) shows that 
the damping effect of the phase fluctuation on the entan- 
glement recently reported when there is no driving ra- 
diation could have been compensated by decreasing the 
amplitude of the driving radiation in case Q 2 < 2Tj and 
increasing when Vt 2 > 2T-f. For arbitrary amplitude of 
the external driving radiation, the effect of phase fluctu- 
ation turns out to be significant when the atoms initially 
have equal probability to be either in the lower or the 
upper energy level. Earlier communications indicate that 
strong entanglement is observe in the vicinity of r\ = 0.1, 
that is when initially the atoms have nearly 48% to be 
in the upper energy level. Putting these together entails 
that the partial initial preparation may cause the appear- 
ance of entanglement in this case since it is not possible 
to maintain the atoms at 50% probability to be in ei- 
ther state even if ones wishes to do so. This implies that 
phase fluctuation may indirectly be a source for the man- 
ifestation of quantum features. It is worth noting that a 
similar argument can be applicable to demonstrate the 



increment of the mean photon number with phase fluctu- 
ation, contrary to the expectation and recent report for 
arbitrary values of 77 and when £1 = [|[ . 

It is also tempting to consider what will happen if it 
is possible to maintain a definite phase between the up- 
per and lower energy levels during the preparation rather 
than phase fluctuations we have assumed in this work. 
The assumption that one has a complete control over the 
phase preparation via locking the laser deployed for ini- 
tial prepapration to a particular known phase leads to 
entirely a different picture. Following the same approach 
without taking an average over the phase shows that both 
the mean photon number and quantum features of the 
generated radiation exhibit oscillatory nature. From the 
face value, this evokes that if we can properly harness 
this mechanism, it may lead to enhancement in the de- 
gree of nonclassicality of the radiation. To see the actual 
effect in this regard, a detailed complete analysis is re- 
quired. Therefore, based on this study, it can be observed 
that the inclusion of the phase through initial prepara- 
tion of the atoms can significantly enrich the potential of 
the nondegenerate three-level cascade laser as a source 
of nonclassical light which requires an indepth further 
analysis. 
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